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Magnetic impurities in three-dimensional Dirac and Weyl systems are shown to exhibit a fascinat¬ 
ingly diverse range of Kondo physics, with distinctive experimental spectroscopic signatures. When 
the Fermi level is precisely at the Dirac point, Dirac semimetals are in fact unlikely candidates 
for a Kondo effect due to the pseudogapped density of states. However, the influence of a nearby 
quantum critical point leads to the unconventional evolution of Kondo physics for even tiny devi¬ 
ations in the chemical potential. Separating the degenerate Dirac nodes produces a Weyl phase: 
time-reversal symmetry-breaking precludes Kondo due to an effective impurity magnetic field, but 
different Kondo variants are accessible in time-reversal invariant Weyl systems. 


The last decade has seen enormous activity in the re¬ 
search of electronic topological states of matter, sparked 
by the discovery of topological insulators (TIs), protected 
by time-reversal symmetry (TRS) [THS]. More recently 
it was appreciated that gapless bulk states can also be 
topologically nontrivial leHU]. Experimental candidates 
for such three-dimensional (3D) Dirac semimetals (SMs), 
such as NasBi and Cd 3 As 2 , are now emerging [T5H2Q]. 

Dirac SMs possess degenerate Dirac nodes at isolated 
points of the Brillouin zone, around which electronic ex¬ 
citations behave as 3D massless Dirac fermions. As such, 
Dirac SMs are 3D counterparts of graphene m, whose 
doubly-degenerate Dirac cones are protected by TRS and 
lattice inversion symmetry (IS). 

The degeneracy can be lifted by perturbations break¬ 
ing either TRS or IS. Nodes of opposite chirality are then 
separated in momentum space or energy in the resulting 
topological Weyl phase. This separation produces re¬ 
markable stability to interactions m, and unique phys¬ 
ical properties such as the ‘chiral anomaly ’ m and sur¬ 
face Fermi arcs [9]. Several Weyl SM materials have re¬ 
cently been discovered [22112^ - including TR-invariant 
systems such as TaAs and NbAs [25ll27] . 

In this paper we discuss Dirac and Weyl systems 
with isolated magnetie impurities (e.g. transition metal 
adatoms), the scattering from which producing a diverse 
range of Kondo physics whose characteristic experimen¬ 
tal signatures fingerprint the topological host. We show 
that the host density of states (DoS), and the symmetries 
broken by perturbations, play a key role in determining 
the kind of Kondo effect that can take place [28] . Pertur¬ 
bations can lead to Dirac metals and insulators or vari¬ 
ous Weyl variants; in each case, different Kondo physics 
manifest, ranging from standard metallic, through quan¬ 
tum critical, to pseudogap Kondo. Fig. summarizes 
these findings. In all Kondo phases, we uncover an un¬ 
conventional evolution of the Kondo temperature when 
the Fermi level is tuned in the vicinity of the Dirac point 
— see Fig. 

The theoretical framework we develop applies to im¬ 
purities in Dirac materials such as NasBi and Cd 3 As 2 ; 
we show that strong particle-hole asymmetry (which can 
be on either the impurity or in the high-energy part of 
the host DoS) can stabilize a pseudogapped Kondo state. 
On doping the Fermi level away from the Dirac point, the 
Kondo effect ubiquitously appears, although influence of 


a nearby quantum critical point is important. In Weyl 
systems, we show that the specific way in which the nodes 
are separated in momentum space is important. When 
it is achieved by TRS breaking (e.g. due to magnetic or¬ 
dering, as proposed in pyrochlore iridates 0 ), the Kondo 
effect is suppressed. When the nodes are separated due to 
IS breaking (e.g. in the newly-discovered TRS-invariant 
but non-centrosymmetric Weyl materials TaAs and NbAs 
[25ll27] i. pseudogap Kondo can again be realized. Any 
given microscopic Dirac/Weyl realization can be classi¬ 
fied according to its symmetries, and can be compared 
to one of the generic situations discussed in this Letter. 
Alternatively, details of the microscopic structure can be 
inferred from the type of Kondo effect observed. 

Bulk model and symmetries: The clean, unperturbed 
Dirac SM is described by the minimal Bloch Hamiltonian, 

TLbO^) = vf Tz a , (1) 

where a and f are Pauli matrices acting respectively in 
spin and orbital space, and r’i? is the effective Fermi ve¬ 
locity. The pristine Dirac SM possesses both TRS and 
IS — meaning T'HbO^)T~^ = TLb{—^) in terms of the 
time-reversal operator T = tq 0 {icFy)K (with complex 
conjugation iF); and = 7iD(—k), with in¬ 

version operator V = 

Bulk perturbations: The description of real systems ne¬ 
cessitates inclusion of perturbations to Eq. 0 ; in partic¬ 
ular, the node-separation in Weyl systems is by defini¬ 
tion a result of these perturbations, and therefore they 
must be considered explicitly. To leading order, pertur¬ 
bations are of form 6H = • r + aofo^ ^ (b ' a P • 

Microscopic models of Dirac/Weyl systems map onto 
TLbO^) + low energies. Depending on the param¬ 

eters a, ao, b, and bo, the perturbations can leave the 
Dirac system TRS and IS invariant, break either TRS or 
IS to produce a Weyl system, or break both TRS and IS 
[29] . Representative examples are discussed below; see 
also Fig. 

Impurity model: We model the magnetic impu¬ 
rity as a single correlated quantum level, coupled lo¬ 
cally in real-space to conduction electrons of the host 
Dirac or Weyl system. The Andersonian Hamiltonian 
|3Q| is H — Hhost T -^imp T -^hyb: whcrC -f^host — 

/ (Cys ^^(k)['^D(k) + (5^]T(k) in terms of conduction 
electron operators T(k) in r- and a-space. The im- 
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FIG. 1. Summary of possible perturbations to the Dirac semimetal, SHi to SHq (top to bottom), and their effect on the 
dispersion (sketched in the space of kx, ky and E for fixed kz), spin-resolved density of states pcr{(^), and scattering T-matrix 
spectrum t(uj) at T = 0 (left to right). Shown for band cutoff D = vf = ^ and magnetic impurity interaction strength U = 0.3. 
Impurity asymmetry 77 = —2elU = 1 and impurity-host coupling gr = g = 0.2 for ( 5 i 7 i, 2 , 3 , 4 , 6 , while for H 5 we use 77 = | 
and g = 0.42,0.479,0.52 for the ALM, critical ACR, and ASC results, respectively. Perturbation strength /i = — | for SHi, 
m = 0.1 for SH 2 , bz = 0.1 for (5 7 / 3 , 4 , 5 , and a = | for SHq. The T-matrix results in the fourth column are calculated for the 
impurity-coupled system using the DoS shown in the third column, obtained exactly by NRG. 
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purity is described by i^imp = 

with da- a spin a =t, i operator for the impurity level. 
T] = — 2e/[/, characterizes the impurity particle-hole 
(PH) asymmetry. The impurity-host coupling is speci¬ 
fied by i^hyb = jZa,r = 0) +H.C.. The CXact 

retarded impurity Green function, after integrating out 
bulk fermions, reads = [(2; —e)t —A(2;) — 

where the hybridization matrix is given by Aa^a'i^) = 

in terms of local host Green 
functions at the impurity position x = 0 , given by 
ghost ^ j Depending 

on the microscopic details of how the impurity is em¬ 
bedded in the Dirac/Weyl host material, the coupling 
gar iiia.y depend explicitly on a and r. However, the 
hybridization is a 2 x 2 matrix in spin space, and in¬ 
volves a trace over the pseudospin degree of freedom r. 
Ultimately, the Kondo physics is controlled by the sym¬ 
metries of A{z). Although for concreteness we now set 
gar = g, real systems should be comparable to one of the 
distinct cases discussed here. 

I](z) is the interaction self-energy, which contains non¬ 
trivial correlation effects due to the impurity. We ob¬ 
tain it exactly from the numerical renormalization group 
(NRG) [3T1|32], whose input is A(^). 

T-matrix: The T-matrix describes scattering in the 
Dirac or Weyl system due to the impurity, 0^°®^(^) = 
ghost ^ with elements Tar.a'r'iz) = 

g'^g^,G^^^,{z). Related to the impurity DoS, the T- 
matrix can be probed locally by scanning tunneling spec¬ 
troscopy. Quasiparticle interference is also sensitive to 
the T-matrix EH El, as is resistivity [35]. Here, we con¬ 
sider its spectrum, t{uj) = — Tr T{uj ^ i0+). 

Kondo physics in Dirac systems: The energy eigenval¬ 
ues of are two-fold degenerate. The local host 

DoS in 3D, p^{uj) = -UmEr + z0+) oc is 

therefore quadratic at low energies, since the Dirac point 
is at the Fermi level. Furthermore, A{z) is diagonal in 
spin-space, and so the effective impurity problem falls 
into the pseudogapped Kondo class [36f[38] . with intact 
TRS. With 7 ^ = 1 on the impurity, the entire system is 
PH symmetric; the depleted DoS strictly precludes the 
possibility of a Kondo effect. Below temperature/energy 
scales ^ U, a spin-^ local moment (LM) forms which 
remains unscreened down to T = 0. 

When PH symmetry is broken, 77 7 ^ 1, the impurity can 
be Kondo screened — there is a quantum phase transition 
from the asymmetric LM (ALM) phase to an asymmetric 
strong coupling (ASG) Kondo phase, on increasing the 
impurity-host coupling with a nontrivial interacting 
critical point (AGR) at gc [37]. The Kondo temperature 
in the ASG phase vanishes as Tk ^ Id — 9c\ on approach¬ 
ing the critical point g ^ g^. Incipient RG flow in the 
vicinity of AGR arises in both ASG or ALM phases for 
small \g — gc\ ^ U. We have confirmed this scenario 
for the Dirac SM using NRG, albeit that a rather large 
coupling g is required to access ASG. 

The effects of perturbing the pristine Dirac SM are 
discussed below in relation to Fig. (top to bottom for 


different perturbation types), showing in particular the 
exact T-matrix in the right column. In Fig. [^we study 
the unconventional evolution of the Kondo temperature 
Tx, in all cases where the Kondo effect manifests. 

(I) dHi = —jiTo G (Jo corresponds to a chemical po¬ 
tential and does not break TRS or IS — see first row 
of Fig. The low-energy DoS Pa{^) oc (cj -h /i)^ is fi¬ 
nite at the Fermi level, implying that the impurity spin 
is screened on the lowest energy scales. For large /i, we 
indeed find a finite T^, which scales in the ‘standard’ 
way for magnetic impurities in metals m, 

Tif ~ ; J oc g^[U7]{2 - T])]-^ , (2) 

with J the effective Kondo exchange. The top-right panel 
of Fig.l^shows the spectrum of the T-matrix due to impu¬ 
rity scattering at T = 0, exhibiting the classic three-peak 
structure. The central Kondo resonance (of width pro¬ 
portional to Tk) embodies enhanced spin-flip scattering 
at low energies, and is pinned at the Fermi level ccljp? 
by the Friedel sum rule [30|. Note the pronounced PH 
asymmetry of the spectrum, due to the asymmetric DoS. 

A metallic Kondo effect would imply the scaling 
InTx ^ meaning the Kondo temperature be¬ 

comes very small (and effectively unobservable in experi¬ 
ment) when the Fermi level is near the Dirac point. How¬ 
ever, this analysis breaks down at small p due to the 
nearby AGR quantum critical point. 
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FIG. 2. Evolution of the Kondo temperature in the three sit¬ 
uations where the Kondo effect can manifest. Tk vs chemical 
potential p shown for the Dirac metal with 6 Hi, and the line- 
node Weyl semimetal with 6 Hi -h 6 H 5 , in the top and middle 
panels, respectively. Bottom panel shows Tk vs a in the Weyl 
metal with SHq. Dotted lines show the powerlaw crossovers. 
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For an asymmetric impurity 7 ^ < 1, the impurity is 
Kondo-screened in the ASC phase for g > Qc^ even at 
/i = 0. The Kondo temperature in fact evolves smoothly 
from this finite /i = 0 value to Eq. [^on increasing /i > 0. 
This is shown in the top-right panel of Fig.(solid lines). 
By contrast, Tk = 0 in the ALM phase for g < Qc c^t 
/i = 0 , implying Tk ^ 0 as /i ^ 0 (see dashed lines). 
The infiuence of the critical point therefore extends well 
into the ‘standard’ Kondo phase. Indeed, the finite-/i > 0 
Kondo phase is separated into a ‘high-T^ regime’ and a 
‘low-T/^ regime’ along Tk ^ \g — gc\ (dotted lines, Fig.[^, 
depending on the parent /i = 0 phase. 

PH symmetry is broken by both f] I and /i 7 ^ 0, 
which can either reinforce or partially cancel. Varying 
/i therefore changes the effective proximity to the critical 
point. For a given 77, this leads to a pronounced asymme¬ 
try between g > 0 and /i < 0 , as seen in the top panels in 
Fig. Indeed, increasing \g\ can drastically reduce Tk 
by tipping the system from the high-T^ (ASC parent) 
regime to the low-T^ (ALM parent) regime. Only for 
very large |/i| is Eq. [^recovered. 

(II) 5 H 2 = mix C) do opens a gap — see Fig. (sec¬ 
ond row). There is no Kondo screening at PH symmetry 
because there are no low-energy electronic degrees of free¬ 
dom to build up the Kondo singlet. The ground state is 
therefore a free LM. Although the T-matrix for this sys¬ 
tem at T = 0 contains Hubbard satellites at energies e 
and ( 1 / + e) due to charge ffuctuations, it is fully gapped 
at low energies, as verified by NRG. For PH asymmetric 
systems, it is known from semiconductors/insulators [39] 
that a first-order transition between LM and Kondo sin¬ 
glet states arises on increasing g. In principle, mid-gap 
poles in the T-matrix could appear. But in the present 
context of impurities in Dirac insulators, this scenario 
seems unlikely, due to the combination of strong PH 
asymmetry and large g which we find are required. 

TRS-broken Weyl systems: 

(HI) SHs = To 0 b ' a splits the two degenerate Dirac 
theories in momentum space — see Fig. (third row). 
Without loss of generality, we choose b = bz^ such that 
gimp diagonal. Although the low-energy DoS for a 
and I has the same quadratic form to leading order, the 
bulk TRS-breaking leads to an effective magnetic field 
on the impurity, A>[^>[^(0) — A||(0) oc bz. The impurity 
is fully polarized on the temperature/energy scale of bz, 
and no Kondo effect can occur. This is manifest in the 
(spin-summed) T-matrix as sharp peaks, each of weight 
^ 4^^. At low energies, the spectrum must decay with 
the same power as the hybridization, t{uj) ^ 

(IV) 6H4 = Tx 0 b ' a: one of the two Dirac theories is 
gapped out, while the other theory acquires a line node 
— see Fig. [^(fourth row). This line-node Weyl SM has a 
linear DoS gOj. TRS-breaking induces a finite effective 
impurity field, and Kondo is again quenched. The low- 
energy behavior of the T-matrix is now linear, t{uj) ^ jcjl. 

TRS-invariant Weyl systems: 

iy) 6H^ = Ty 0 b ' a realizes a line-node Weyl SM as 
for case (IV), but without an impurity field, A>^^(z) = 
A 4 ,^( 2 ;). The low-energy DoS is linear [40], pai^) oc |c(;|. 


A magnetic impurity in this line-node SM is therefore a 
good candidate to realize pseudogap Kondo physics: it 
is more likely accessible than in the Dirac SM case be¬ 
cause the DoS vanishes less rapidly near the Fermi level. 
The ASC phase can therefore be reached with moder¬ 
ate impurity PH asymmetry 77 and smaller impurity-host 
coupling g. The fifth row of Fig. contains NRG results 
for the T-matrix in the ALM/ASC phases, and at the 
ACR critical point, each showing distinctive signatures. 

As with case (I), the pseudogap Kondo critical point 
has a significant infiuence on the ‘regular’ Kondo phase at 
finite /i, demonstrated by the unconventional evolution of 
Tk in Fig. [^(middle panels). For 77 < 1 and /i > 0, there 
is again a separatrix between high- and Iow-Tk regimes 
following Tk ^ |/i|. At large g, the Kondo temperature 
varies more slowly than in the Dirac SM, InT^ r\j -i/H 
However, the competition between asymmetry coming 
from the impurity (77 7 ^ 1 ) and host {g 7 ^ 0 ) is more 
finely balanced due to the linear DoS. This is starkly ev¬ 
ident for g < 0 , where a nontrivial intermediate regime 
of two-stage Kondo screening is found between Tk ^ \g\ 
and Tk ^ |/7|^ (dotted lines), with ^ 2.8 extracted 

from NRG. In this regard, the TRS-invariant line-node 
SM is similar to graphene, where equivalent impurity ef¬ 
fects have previously been found m- 

(VI) 6Hq = aiz 0 do splits the two degenerate Dirac 
theories in energy — see sixth row. Fig. The DoS is 
nonvanishing everywhere, and one expects regular Kondo 
physics (see e.g. the spectrum of the T-matrix, plotted 
for large perturbation strength a = |). For small a, we 
again see the infiuence of the parent a = 0 pseudogap 
Kondo quantum critical point, as shown in Fig. (bot¬ 
tom panel). The situation here is distinct from cases (I) 
or (V), since the perturbation SHq does not introduce 
additional PH breaking. We find from NRG a new and 
nontrivial powerlaw separatrix between high- and low- 
Tk regimes, Tk ^ Ini'", with ( = 1.8. The underlying 
physical explanation for these power laws is left for fu¬ 
ture study. 

Conclusion: Dirac and Weyl systems are shown to con¬ 
stitute a remarkably rich playground for Kondo physics, 
potentially exhibiting essentially the full range of single¬ 
impurity effects. In particular, the unconventional evo¬ 
lution of the Kondo temperature on varying chemical 
potential, and the unusual role of PH asymmetry, are 
‘smoking gun’ experimental signatures. 

We anticipate further interesting physics in topologi¬ 
cal Kondo lattice systems, and those with several mag¬ 
netic impurities where pseudogap Kondo, RKKY, two- 
impurity Kondo, and disorder compete. 
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